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In this paper, we apply the first-principle theory to explore how the electronic structures of 
armchair graphene nanoribbons (AGNRs) are affected by chemical modifications. The edge addends 
include H, F, N, NH2, and NO2. Our theoretical results show that the energy gaps are highly tunable 
by controlling the widths of AGNRs and addends. The most interesting finding is that N-passivated 
AGNRs with various widths are metallic due to the unique electronic features of N-N bonds. This 
property change of AGNRs (from semiconducting to metallic) is important in developing graphene- 
based devices. 

PACS numbers: 73.22.-f, 73.20.Hb, 72.80.Rj, 73.63.Bd 



I. INTRODUCTION 

Graphene, a single atomic layer of graphite with a hon- 
eycomb crystal structure, has attracted a great number 
of research activities^ Its structural properties, electri- 
cal conductance and quantum Hall effects have been in- 
vestigated for making novel nanoelectronic devices 
Due to the linear energy dispersion relation near the 
Dirac points, graphene is an interesting conductor in 
which electrons move like massless Dirac- fermions.^ It 
is now possible to make graphene nanoribbons (GNRs) 
with various experimental methods such as tailoring via 
a scanning tunneling microscopy tip,^ exfoliating from 
high oriented pyrolytic graphite, or graphitizating 
SiC wafers.^ The energy gaps of GNRs with several tens 
of nanometers can be measured because the growth of the 
energy gap is inversely proportional to the GNRs' width. ^ 
There are two typical types of GNRs according to their 
edge configurations, either armchair or zigzag. Zigzag 
edge GNRs are metallic without considering the free- 
dom of spin because the two edge states are degenerated 
at the Fermi level. Hydrogen-saturated armchair GNRs 
(AGNRs), on the other hand, are semiconductors i^ii^ 
The electronic transport and magnetic properties of 
zigzag GNRs have also been studied by several research 

groups i^iiiii^ii^ 

Because of the high ratio between edge and inner 
atoms, the electronic structure of narrow GNRs are sen- 
sitive to various edge addends. Individual carbon atoms 
on graphene edges are only bounded to two neighboring 
carbon atoms and a dangling carbon bond offers a re- 
markable opportunity for altering GNR's electronic prop- 
erties. This edge modification can be implemented by 
attaching various atoms or functional groups to these 
dangling carbon atoms. To the best of our knowl- 
edge, edge modified AGNRs, except for the hydrogen- 
passivated GNRs, have not been systematically exam- 
ined using Density Functional Theory (DFT). In this 



paper, the electronic structures of AGNRs, with vari- 
ous addends including H, F, N, NO2 and NH2, are thor- 
oughly investigated based on first-principle calculations. 
Our results show that AGNRs can be either semicon- 
ducting or metallic by changing edge chemical addends. 
This remarkable characteristic is very useful in making 
graphene-based molecule computing or sensing devices. 
The detailed analysis of the electronic structures can also 
be used to understand the underlying microscopic mech- 
anisms of these edge effects. 



II. MODEL AND METHOD 

The electronic structure calculations are carried out by 
employing the Vienna ah initio simulation package with 
local density approximation. ^^^^^ The electron- ion inter- 
actions are described based on the projected augmented 
wave (PAW) and the frozen core approximation.^^ The 
energy cutoff is set to be 400 eV. Following the previ- 
ous definition of AGNRs, ^ we choose AGNRs with width 
W=18, 19, and 20 (corresponding to 3p, 3p+l, and 
3p+2, respectively. Here, p is an integer) as examples to 
study three kinds of AGNRs. The schematic of a sample 
AGNR with width of W=18 is shown in Figure 1. Vari- 
ous addends, such as H, F, and N atoms, are connected 
to individual edge carbon atoms. In the NH2 and NO2 
passivated cases, each edge carbon atom bonds to one 
NH2 or NO2 group alternatively (its neighboring carbon 
atom is saturated by an H atom). In our calculations, 
all atomic positions are allowed to relax. The conver- 
gence tolerance in energy and force is 2xl0~^ eV and 
0.02 eV/A, respectively. The vacuum layers between two 
neighboring AGNR sheets and AGNR edges are set to be 
12 A and 15 A wide, respectively. In a typical calculation, 
a one-dimensional periodic boundary condition along the 
edge direction is imposed. We choose the k-point sam- 
pling consisting of 15 uniform k-points together with the 




FIG. 1: (Color online) Schematic of a sample AGNR with 
width W=18. The black dots correspond to carbon atoms 
while the blue ones represent the addends (labeled with X, 
where X can be H, F, N, 2H, 2F, NO2, and NH2 molecules). 
The distance between two dash-dotted lines is the lattice con- 
stant of one-dimensional unit-cell in the DFT calculations. 



TABLE I: The calculated energy gaps (in eV) of X-AGNRs 
with width w=18, 19, and 20. Here, "M" stands for the metal- 
lic case without Egap- 
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^Because there is not much difference in bandgaps when nar- 
row graphenes are selected, only the formation energy of AGNRs 
with width w=18 is considered. The formation energy is defined 
as, AEf = ^{Etot - Etare " E^xMx), where Etot, Etare, n, 
nx, and fix is the total energy of system, the energy of un- 
passivated 18-AGNR, the number of chemical groups, the num- 
ber of atoms of addends X and the chemical potential of X, re- 
spectively. We choose the chemical potential of H, F, N, and O 
as the binding cohensive energy per atom in H2, F2, N2, and 
O2 molecules. This definition would be a appropriate mesure- 
ment of the stability in passivated AGNRs. The formation energy 
of NO2 and NH2 passivated 18-AGNRs per group is defined as 
AEf = \ {Etot-Ebare - E ^xMx -2 X AEf{H)), where AEf(H)) 
is the formation energy of H-passivated 18-AGNR. 

r point. 



III. RESULTS AND DISCUSSIONS 

We first choose the well-studied hydrogen passivated 
AGNRs^'^^ as a basis to validate our numerical results. 
The calculated band structures of three kinds of hydro- 
gen passivated AGNRs are plotted in Figure 2 (a). The 
corresponding bandgaps together with the optimized C- 
C distances at the edge are listed in Table I. Our re- 
sults show that all three kinds of AGNRs are semicon- 
ductors with direct band gaps at the F point. The figure 
clearly shows that the AGNR with width W = 3p ^ 1 
has the largest energy gap, while the AGNR with width 
W=3p+2 has the smallest one. Our calculated bandgaps 
are 0.41, 0.61, and 0.10 eV for the H passivated AGNRs 
with width W=18, 19 and 20, respectively. These results 
agree well with the previously reported result^^. 

Next, we examine the electronic structures of AG- 
NRs with F atoms passivated at both edges. The cal- 
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FIG. 2: (Golor online) (a) and (b) show the band structures 
of H and F-passivated AGNRs with width W= 18, 19, and 
20, respectively, (c) The projected density of states(PDOS) of 
H (upper panel) and F (lower panel) passivated AGNRs with 
width W=18. 



culated band structures are plotted in Fig. 2 (b) and 
the bandgaps are listed in Table T Similar to the H- 
passivation cases, all F-passivated AGNRs are semicon- 
ductors. However, the effects of GNR width on bandgaps 
is clearly different. For instance, the F-passivation makes 
the band gap of the AGNR with width W=18 (3p) to 
be the largest (0.61 eV), while the gap of the AGRN 
with width W=20 (3p+2) is the smallest (0.14 eV). It 
is important to note that the bandgaps of the AGNRs 
(width W=20) with the smallest gaps remain as the 
smallest even after passivation. However, AGNRs with 
the largest bandgaps behave differently and depend on el- 
ements used for passivation. We also would like to know 
what happens if individual edge carbon atoms connect 
to two H or F atoms instead of one H or F atom. In- 
tuitively, the sp^ hybridization of edge carbon atoms is 
converted into sp^ hybridization. The effective AGNR 
width is therefore reduced by 2. In other words, the 
effective width becomes 17 for both double H- and F- 
passivated AGNRs when their initial width is W=19. 
Because the effective width of 17 fits 3p+2, the corre- 
sponding double H- and F-passivated AGNRs have the 
smallest energy gap (about 0.1 eV). The largest gaps are 
0.70 eV for double H-passivated AGNRs with W=17 or 
effective W=16 (3p+l). Similarly, 0.63 eV for of double 
F-passivated AGNRs with W=20 or effective W'=18 (a 
3p case). These results are consistent with the width de- 
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pendent energy gap features in single F- or H-passivated 
AGNRs. 

It is instructive to compare the DFT results with the 
tight-binding approximation (TBA) data. Wang et al. 
have reported that energy gaps can be tuned by changing 
edge hopping parameters.^'' Son et al have also showed 
that the energy gaps of H-passivated AGNRs can be 
reproduced via TBA simulations. In their simulations, 
edge hopping parameters can be increased by decreasing 
the bond length of C atoms at edges. ^ Although the ge- 
ometric deformations of carbon networks are almost the 
same between H-passivated GNRs and F-passivated ones, 
the bandgaps surprisingly show distinct dependence on 
the width of AGNRs. According to the TBA parameters 
adopted in these analytical calculations of the edge mod- 
ified GRNs,^^ the iVpp hopping integrals, representing the 
coupling between orbital of two neighboring carbon 
atoms, increases as the C-C distance decreases. However, 
this statement is invalid when the bond length substan- 
tially decreases so that the hybridization between two 
carbon atoms becomes perturbed. From the PDOS of 
the H-passivated AGNRs, we can see that the peaks of pz 
orbital of edge C atoms almost coincide with that of the 
inner C atoms below the Fermi level [see Fig2(c)]. The 
TTpp bonding remains unchanged during the H-passivation 
of edge C atoms. The modification of energy gaps caused 
by H atoms could be explained based on the geometry 
deformation. Interestingly, the Pz orbital of both edge C 
atoms and F addends have peaks at —8.5 eV below the 
Fermi energy. The peak position of the edge C atoms, 
however, decreases significantly compared to the inner 
C atoms from —6.0 to —2.3 eV below the Fermi energy, 
which is still within the energy range of tt bonding in F- 
passivated AGNRs. That is to say, the Pz orbital of edge 
C atoms in the F-passivated cases form chemical bonds 
with the neighboring C atoms and the F atoms. Conse- 
quently, the hopping parameters decrease. If we choose 
the hopping parameter to be 84% of -2.7 eV ( the typ- 
ical hopping integral of sp^ carbon system), the energy 
gaps of AGNRs with width W=18, 19, and 20 are 0.65, 
0.35, and 0.13 eV, respectively, which agree excellently 
with results obtained by DFT. This observation demon- 
strates that the electronic structures of AGNRs can be 
controlled by edge chemical modifications. 

NO2 and NH2 radicals are widely used as anchor- 
ing or side groups for making single molecular junctions 
in molecular electronics. ^^'^^^^^ In the next calculations, 
NO2 and NH2 are attached to AGNRs and we investigate 
how these addends impact GNR electronic structures. 
Figs. 3 (a) and (b) show their band structures. Figs. 3 
(c) and (d) show their optimal geometric structures. The 
NH2 group lies in the AGNR plane for NH2 passivated 
GNRs, while, for NO2 passivated GNRs, the dihedral an- 
gle is 59° between AGNR and NO2 plane. According to 
the projected DOS in Figure 4, their electronic structures 
also exhibit different features. For the NH2-passivated 
AGNRs, the nitrogen atom forms one sp^ orbital, which 
bonds to two hydrogen and one edge carbon atom. The 




FIG. 3: (Color online) Band structures of NO2-, and NH2- 
passivated AGNRs with width W = 18, 19, and 20 are show 
in (a) and (b), respectively. The edge relaxed geometric struc- 
tures of NO2 and NH2 are shown in (e) and (f). The relaxed 
bond lengths from ai to ae are 1.36, 1.40, 1.39, 1.37, 1.42, 
and 1.37 A, respectively, and the bond angles from 71 to 74 
are 122.8°, 120.6°, 117.5°, and 126.6°, respectively 



remaining p^ orbital of the nitrogen atom bonds nicely 
to the Pz orbital of the edge carbon. For the NO2 case, 
the TT orbital in the NO2 plane does not bond well to the 
sp^ framework of AGNRs because the NO 2 plane twists 
out of the AGNR sheet. Comparing to the bandgaps 
of H- and F-passivated AGNRs, it is interesting to note 
that the bandgap of NO2 is close to that of H, while the 
bandgap of NH2 is similar to that of F. The fundamen- 
tal cause is attributed to the infiuences of the carbon p^ 
orbital. Neither H nor NO2 forms a bond with the car- 
bon Pz orbital, but both F and NH2 facilitate bonding 
directly to carbon pz orbital. Hence, it is expected that 
substantial charge transfer occurs together with the non- 
negligible change in the hopping integral near the edges 
of both F- and NH2 -passivated AGNRs. 

Very recently, N-rich Oligoacences were considered to 
be candidates for n-channel organic semiconductors.^^ N- 
doped carbon nanotubes and nanofibers have also been 
synthesized and investigated intensively over the past 
years ^^^A^iiS^iS^ These investigations motivate us to ex- 
amine the N-passivated AGNRs. According to the calcu- 
lated formation energy in Table 1, it is possible to attach 
one N atom to each C atom along the edges as an ad- 
dend. Interestingly, two neighboring N addends can form 
an N-N dimer and form an optimal structure. The N-N 
bond length is 1.30 A, which is 0.20 A longer than that 
of N-N bond length in N2 molecule. More importantly, 
a semiconductor-to-metal transition is observed for all 
N-passivated AGNRs and this transition is independent 
of AGNR widths. To explore the cause at the micro- 
scopic scale, the band structures of N-passivated AGNRs 
with various graphene widths and partial projected DOS 
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FIG. 4: (Color online) (a) Partial density of states (PDOS) 
of the edge C, and N in NHs-passivated AGNRs. (b) PDOS 
of the C, N, and O atoms in NO2-AGNRS. Here, W=18. 



on N atoms of AGNR with width W=18 are presented 
in Figs. 5(a) and 5(b), respectively. Clearly, the DOS 
near the Fermi level is impacted by N addends. The 
edge N atom has 5 valence electrons. Two of them stay 
as a lone-pair, which locates in the AGNR' plane with 
the direction pointing outwards. The occupied band la- 
beled with "A" in Fig. 5(a) has the same features as 
those of lone-pair electrons. The other three electrons 
exhibit sp^-like hybridization. Two of the sp^ hybridized 
orbital pair with the neighboring C and N atoms to form 
two a bonds, while the last electron acts as a tt elec- 
tron. We would like to emphasize that the metallic state 
of N-passivated AGNRs is predominantly contributed by 
the edge states of N addends. There are two reasons 
attributed to AGNR's semiconducting to metallic tran- 
sition, i) the N-N bond is mainly composed of 2s and 
2py orbit als. The energy of the band with N 2pz is above 
the Fermi surface, which means that N 2pz electrons fill 
the carbon tt* bands and push the Fermi surface inside 
the carbon's tt* bands, ii) the lone-pair electrons of N 
atoms are mainly in N 2px bands and the top part of 2px 
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FIG. 5: (Color online) (a) Band structures and DOS of N- 
passivated AGNRs with width W = 18, 19, and 20, respec- 
tively. Here, the total DOS and projected DOS on N atoms 
are labeled with solid and dashed lines, (b) Partial projected 
DOS on N atoms of N-passivated AGNRs with width W=18. 
The relaxed geometry is shown in the inset in Fig.(b) 



bands pass the Fermi surface. This feature contributes 
to the semiconducting to the metallic transition. Clearly, 
edge N-N addends extend the sp^ network of AGNRs and 
thus the effective width (or the number of sp^ dimer lines) 
of N-passivated AGNRs increases by 2 (coming from two 
edge sides). Apart from "A" band, the band structures of 
N-passivated AGNRs with width W are similar to those 
of H or F-passivated AGNRs with width W -\- 2 except 
that the Fermi level shifts upwards and cross the lowest 
unoccupied bands, which leads to the semiconductor-to- 
metal transition. 



IV. CONCLUSION 

In summary, we perform the first-principle calculations 
on the electronic structures of AGNRs with various edge 
chemical modifications. We find that the energy gap can 
be tuned to be either metallic or semiconducting by at- 
taching different chemical functional addends. After AG- 
NRs are passivated by H, F, NH2, and NO2, the GNRs 
with the smallest bandgaps (i.e. W=20) remain to have 
the smallest bandgaps. However, AGNRs that have the 
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largest bandgaps behave differently and the gaps depend 
on the type of elements used for passivation. H and F 
have similar bandgaps to NO2 and NH2, respectively, 
which can be explained based on the bonding of the car- 
bon pz orbital. The fascinating observation is that AG- 
NRs passivated by attaching N atoms can change from 
semiconducting to metallic material regardless of AGNR 
width. This result is analyzed in detail based on its 
electronic structure. The unique electronic structure of 
the N-N bond, which is responsible for semiconductor-to- 
metal transition, is explained based on its 2s and 2py or- 
bital, and the lone-pair of 2px. This tunable character of 
AGNRs is useful in developing graphene-based molecule 
electronic devices in the near future. 
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